BUTP-98/1 



The index theorem in QCD with a finite cut-off Q 



Peter Hasenfratz, Victor Laliena, 
Institute for Theoretical Physics 
University of Bern 
Sidlerstrasse 5, CH-3012 Bern, Switzerland 



Ferenc Niedermayeif] 
Paul Scherrer Institute, CH-5232 Villigen PSI, Switzerland 



February 2008 



Abstract 

The fixed point Dirac operator on the lattice has exact chiral zero modes 
on topologically non-trivial gauge field configurations independently whether 
these configurations are smooth, or coarse. The relation riL^TiR — Q FF , where 
til (nit) is the number of left (right)-handed zero modes and Q FP is the fixed 
point topological charge holds not only in the continuum limit, but also at 
finite cut-off values. The fixed point action, which is determined by classical 
equations, is local, has no doublers and complies with the no-go theorems 
by being chirally non-symmetric. The index theorem is reproduced exactly, 
nevertheless. In addition, the fixed point Dirac operator has no small real 
eigenvalues except those at zero, i.e. there are no 'exceptional configurations'. 
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In the continuum the Dirac operator of massless fermions in a smooth back- 
ground gauge field with non-zero topological charge has zero eigenvalues. The cor- 
responding eigenfunctions are chiral and the number of left- and right-handed zero 
modes are related to the topological charge of the gauge configuration: n^—n^ = Q 
0| . It is believed by many that this theorem has important consequences on the 
low energy properties of QCD. A possible intuitive picture of a typical gauge con- 
figuration in QCD is that of a gas or liquid of instantons and anti-instantons with 
quantum fluctuations [^|. For such configurations, the Dirac operator is expected to 
have a large number of quasi-zero modes which could be responsible for spontaneous 
chiral symmetry-breaking 

Typical gauge field configurations contributing to the path integral are not 
smooth, however. In addition, the regularization which is necessary to define the 
quantum field theory breaks some other conditions of the index theorem as well. 
Standard lattice formulations, for example, violate the conditions of the index theo- 
rem in all possible ways: the topological charge of coarse gauge configurations is not 
properly defined and the Dirac operator breaks the chiral symmetry in an essential 
way. In the continuum limit these problems disappear, while close to the contin- 
uum some trace of the index theorem can be identified Q. In general, however, the 
expected chiral zero modes are washed away and, even worse, unwanted real modes 
occur as lattice artifacts ('exceptional configurations' f\ ) ||. 

We are going to show here that a lattice formulation of QCD working with the 
fixed point (FP) action (classically perfect action) [6-11] and FP topological charge 
12, [l3| which are defined in the context of Wilson's renormalization group theory 
14], solves all these problems: the index theorem remains valid on the lattice even 
if the cut-off is small (the resolution is poor). We find it amazing that a theorem 
which is based on differential geometry and topology finds its validity in a context 
where none of these notions seem to be defined. 



The FP action in QCD is determined by classical equations. Many aspects of 
these equations and of their solutions have been studied earlier. In particular, it 
can be shown that the symmetries of the continuum theory (infinitesimal transla- 
tion, rotation, chiral transformation, etc.) are realized on the solutions of the FP 
equations of motion although the FP lattice action itself is not manifestly invariant 
jll|| . Actually, the index theorem on the lattice can be demonstrated by using these 
FP equations to connect the lattice problem with that of the continuum, where the 
Atiyah-Singer theorem is valid fLOfl , Here we shall follow another way which allows 
to prove the index theorem on the lattice directly and leads to additional results as 
well. 



We shall first give the essential steps of the argument and then present some 
details in the second part of the paper. 

The fermion part of the FP action can be written as Q 

Sf p (^,V^) = ]T^/w(£/)Vv , (l) 

where h nn > (U) is a specific lattice form of the Dirac operator. In eq. ([!]) the colour 
and Dirac indices are suppressed. We consider one flavour, Nt = 1. The FP Dirac 

3 We call a configuration exceptional if, due to some real eigenvalue close to the origin, the 
quark propagator becomes singular when the bare mass is still distinctly different from its critical 
value. 

4 We use dimensionless quantities everywhere. The dimensions are carried by the lattice unit a. 
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operator is local, has no doublers, satisfies the hermiticity property 

ftt = 75 h l5 , (2) 

and, in complying with the Nielsen-Ninomiya theorem Jl5| , it is not chiral invariant. 
This is realized, however, in a very special way — the chiral symmetry breaking 
part of the fermion propagator is given by 

\{Kl>( U )^}^ = R nn'{U), (3) 

where {, } denotes the anticommutator. The hermitian matrix R nn i(U) is trivial in 
Dirac space and it is local. This chiral symmetry breaking term appears only due to 
the non-chiral-invariant blocking transformation of the fermion fields, and not by 
an ad hoc term added to the action. Due to the locality of R the chiral properties of 
the propagating quark are not affected by this breaking term. The precise form of R 
depends on the block transformation whose FP we are considering. There are block 
transformations which give simply R nn ' = const • 5 nn i . When the averaging extends 
over the hypercube only (as it is the case for all the block transformations considered 
so far), then n~n' extends only over the hypercube as well, but in general R has an 
extension of 0(a). Eq. (|^) with a local R is a highly non-trivial property of the FP 
action. Indeed, is the propagator in the background gauge field U, hence it 

is non-local. For a general case, its chiral symmetry breaking part is expected to be 
non-local as well. For the standard Wilson action and for the currently used other 
improved actions R is non-local, indeed. 

Ginsparg and Wilson observed a long time ago that eq. M) with a local R is 
the mildest way a local action can break chiral symmetry |lq ]. They have also 
demonstrated that the triangle anomaly for U = 1 (free fermions) is correctly 
reproduced in this case. This work remained, however, largely unnoticed since 
no solution was known to the remnant chiral symmetry condition in eq. (|^) for 
QCD. The FP action is a solution as we shall demonstrate in the second part of 
this paper (see also in JlQ]). 

To avoid some singularities at intermediate steps of the calculation we introduce 
a small quark mass which will be sent to zero at the end: 

h, m ' (U) — > hnn' (U) = h nn > (U) + m q S nn i . (4) 

Due to eq. (^) the chiral symmetry breaking part of h satisfies the relation 

75}75 = {h)Rh) nnl +m q 5 nn , . (5) 
This equation will play a basic role in deriving the following results: 



I. The expectation value of the divergence of the singlet axial vector current in 
a background gauge field has the form 

<V M Jj(n)) = -$«76 {h' i Rh^) n + (^Rh) nl5 ^j n ) - 2m q (^ 75 ^) . (6) 
In eq. (||) the expectation value is defined as 

{o§, i>, u)) = ^ry f DijDTpo(i>, v, u) ex P {- J2 i>Ji mn m n } , (7) 
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where Z(U) is given by the analogous integral without the 0(-0, ip, U) factor. In the 
m q — > limit ^|, the last term in eq. (^) is dominated by the zero modes (if there 
are any) 

"L "R 

m lim o m q <^ 7 5^n) =E^ , l > W^' ) W -J2^R*(^R(n) , (8) 

IUq ^ 3=1 «=1 

where ^ (n), tpf}{n) are the normalized A = right- and left-handed eigenfunc- 
tions of the eigenvalue problem of the Dirac operator 

J2 h nn'(U)V n , = A*„. (9) 



II. Summing over n in eq. (|fj) one obtains 

= 2Tr( l5 hR) + 2{n R -n L ) 1 (10) 

where the trace Tr is over colour, Dirac and configuration space. The first term is 
a pseudoscalar gauge invariant functional of the background gauge field. Further, 
it can be shown that it is a FP operator. On the other hand, its density behaves 
for smooth fields U as 



(11) 



up to terms which are of higher order in the vector potential and/or derivatives. It 
follows then 

Tr( 75 W?) = Q FP , (12) 

where Q FP is the FP topological charge defined earlier in the Yang-Mills theory Jl3| . 
As argued in |L2], |l3| , Q FP (if it is used together with the FP action) defines a correct 
topological charge: it associates a unique integer number q to any configuration |^|, 
where q satisfies the inequality: action > 87r 2 |g|, with equality for the instanton 
solutions. Eqs. (0,|l|) give the index theorem. 

III. Eq. (j|) constraints the spectrum {A} of the Dirac operator. For renor- 
malization group transformations ('blocking out of continuum' [l7j| ) which lead to 
Run' — 1/ftf • 5 nn ' , where Kf is a constant parameter entering the block transforma- 
tion, the eigenvalues A lie on a circle whose center is on the real axis in the point 
Kf/2 and has a radius Kf/2 (fig. |l|). In general, for a local R, the eigenvalues lie 
between two circles tangent to the imaginary axis - the the circle described previ- 
ously and a smaller one, as shown on the second picture in fig. ^. This spectrum 
excludes the existence of exceptional configurations. 

We go through the steps now leading to the results above. As usual, the axial 
vector current is constructed from the chiral symmetric part of the action, 

It-sym = ^[ft,75]75 = h - ^{h,j 5 }j 5 . (13) 



5 The m q — > limit will always be taken even if it is not indicated explicitly in the following 
equations. 

6 There exist configurations whose topological charge might change under a small deformation. 
This is the case when a small (size O(a)) topological object falls through the lattice. In the saddle 
point equation, which defines the fixed point action, there occur two degenerate absolute minima 
in this situation, [rj, Il2| . In the very moment when the topological charge changes the number of 
zero modes will also change as required by the index theorem. 
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Figure 1: The spectrum of the fixed point Dirac operator: a) for the special case 
Run' = l//tf • $nn' the eigenvalues lie on a circle, b) for a more general set of block 
transformations the eigenvalues lie between the two circles. 

Performing a local chiral transformation one obtains for the divergence of the singlet 
axial vector current 

V M J*(n) = ^n75 {hsYMip) n + (tphsYu) n TsV'™ , (14) 

where V„/(n) - f(n) - f(n - p). Eqs. (0|1§ give 

-(V3„75 (h^Rh^) n + {^Rh) nl ^ n ) (15) 
-2m q (ip n ~f 5 ip n ) . 

The first term on the r.h.s. is proportional to the equation of motion hip = 0, the 
second comes from the remnant chiral symmetry condition, while the last term is 
due to the explicit symmetry breaking. The contribution of the first term in eq. ( |l5[) 
is zero. Indeed, 

1 f 5 

Z ^J ^ (16) 

and partial integration by %p gives then tr7s = 0. Integrating out the fermions in the 
third term in eq. ( |i"5j ) gives 2m q tY(^c l h~n), where the trace is over colour and Dirac 
space. We insert here the closure relation of eigenvectors of the hermitian matrix 
H = /175 . It is easy to show using eq. d|) that the zero eigenvalue eigenfunctions of 
H (if there are any) are chiral and they are zero eigenvalue eigenfunctions of h as 
well. In the m q — > limit only these eigenfunctions contribute leading to eq. (p|). 

Summing over n in eq. (^) , integrating out the fermions in the first term on the 
r.h.s. and using eq. (^|), one obtains eq. ( |l0| ) immediately. 

Consider now the derivation of eq. (|TT|) . We will show first that 

-2m^ n ^ n ) -,^e^F;„(n)F« p (n) (17) 

up to terms which are of higher order in the vector potential and/or derivatives. 
Eq. (fL7|) can equivalently be written as 



{5 S 



>(0) 5A"Jm) 



= ^ aP 6 ab . (18) 

A=0 



5 



We sketch the derivation only and refer to flq] for the details. Using the notation 

(19) 



m.n f 1 ^ ' 



A=0 

and the definition eq. (0), the l.h.s. of eq. ( |l8| ) can be written as 

((-2m q ^„ 75 ^ n ) j;(m)J, h (0))| A=o . (20) 



If the derivatives in eq. (|T^) act on Z(U), or twice on /i, the contribution can be 
shown to be zero. The current in eq. ( J19| ) is a correct definition for the colour 
vector current: it can be shown that V M J M complies with the Noether theorem. 
The current can also be written as 

m+V , (21) 

i,v 

where T a are the colour generators and Ka is local. 



The matrix element in eq. (20) has been considered in |16j when discussing the 
anomaly (see the considerations in Q between eqs. (29-41)). It can be shown Q 
that in eq. ( ^l| ) satisfies the basic sum rules used in the derivation of the anomaly 
in [^6|. The presence of colour in eq. ( p0| ) gives an extra factor 1/2 • 5 a b in eq(41) in 
leading to eq. (p"7|). In the authors first rewrite —2m q ip n j5ip n in eq. ( pp| ) 
using the operator equation for V M J^(n) and show that only the term coming from 
the remnant chiral symmetry condition (proportional to Rh) contributes. This 



gives the result claimed in eq. ([11]). 

We want to demonstrate now that 

Q(U) = Tr( l5 hR) (22) 

is a FP operator as stated in II.. Further, we have to show that the FP Dirac 
operator satisfies the remnant chiral symmetry condition eq. ([|) . For both problems 
we need the classical equations which determine the FP action. The fixed point lies 
in the /3 ~ I/5 2 — > 00 hyperplane. In this limit the path integral is reduced to 
classical saddle point equations. The gauge field part satisfies |?| 

Sf(V) = mm [Sf(U) + T g (V, U)] , (23) 

where T g (V,U) defines the block transformation relating the block gauge field V 
to a local average of the fine variables U. Denote the minimizing configuration in 
eq. (E3|) by U m i n (V). The FP Dirac operator h is defined by the recursion equation 



KV)n B n' B = K{5 nB n'„ (24) 



y^w(J7 Bl j n ) nB „ (h(U m i n ) + KfbfLU^ (U min )uj(U min )) , w(f7 ro i n )L 



In case h 1 is defined, eq. ( p4[ ) is equivalent to the somewhat simpler recursion 
relation |l9| 

B Kt n,n> (25) 
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Here u>(U), which is trivial in Dirac space, defines the way the block fermion fields 
are constructed from the fine fermion fields, while Kf is a parameter of the block 
transformation. The coarse fields live on a lattice with a lattice unit a' = sa, 
where a is the lattice unit before the transformation, while s is the scale of the 
renormalization group transformation. The factor bf is s * , where is the en- 
gineering dimension of the field (=3/2), while the averaging function satisfies the 
normalization condition ^2 n Lo(U — X) nB n = 1- 

Eq. (|25| ) (or eq. can be solved by iteration. If the scale change s of the 

transformation is (infinitely) large, the fixed point can be reached in a single renor- 
malization group step. In this case, the fine field U m i n lives on an infinitely fine 
lattice, and it is smooth. The Dirac operator h over this field is arbitrarily close to 
the massless continuum Dirac operator. Since lu is trivial in Dirac space, the second 
term on the r.h.s. of eq. ( [25|) is chiral invariant, which leads to the remnant chiral 
symmetry condition in eq. (|3j) with R nn > — 1/nf ■ 8 nn >. If s is finite, one can start 
again on an infinitely fine lattice, but the transformation should be iterated. One 
obtains in every iteration step some contribution to the chiral symmetry breaking 
part of h , and a non-trivial, but local R builds up in this case as claimed in 
eq. (||) before. If h^ 1 is not defined (due to zero modes), one can use eq. (|24|) to 
derive the corresponding remnant symmetry condition for h itself with somewhat 
more algebra. 

By adding to the FP action S FP (U) a small perturbation eO(U), the path in- 
tegral in the saddle point approximation leads to the extra contribution eO'(V) = 
eO(U m i n (V). A FP operator is reproduced by the renormalization group transfor- 
mation up to a trivial rescaling related to the dimension of the operator: 

s - d °0 FP (V) = FP (U min (V)) , (26) 

where s is the scale of the RGT and da is the dimension of the operator O fiofl . 
Since the topological charge is a dimensionless quantity, it satisfies the FP equation 

Q FP (V)=Q FP (U min (V)) ■ (27) 

We want to show that the operator Tr^hR) satisfies eq. (p7|). Using eqs. (p|,p4]) 
and the fact that u> is trivial in Dirac space, one obtains the following recursion 
relation for R 

1 6 2 

R{V)n B n' — 3n B n' H / w(^min)n B n-R(^min)nn'W(J7 m i n ) , , . 

Kf K{ tt> B (28) 

With the help of eqs. (|24|,28) one obtains after some algebra that Q(U) in eq. ( p2f ) 
satisfies eq. (p7|). It is, therefore the FP topological charge studied earlier in spin 
and gauge models (l2[ p^| . 

Finally, let us discuss the spectrum of the FP Dirac operator. Let ip be a 
normalized solution of the eigenvalue equation, hip = Aip. Eq. (|B|) (written in 
the form + h^) = h^Rh) implies ReA = \X\ 2 (tf>, Rip). For the special case 
Run' = 1/ K f ' W the eigenvalues A lie on a circle of radius Kf/2 touching the 
imaginary axis, as shown in fig. [j]. In general, the eigenvalues lie between two 
circles tangent to the imaginary axis - the circle described previously and a smaller 
one as shown in fig. |l|. The smaller circle has a finite radius which follows from 
the locality of R. In particular, one can show that for the block transformation 
used in p|, |2C|] 1/nf < (tp,Rip) < 2/n{ hence the radius of the smaller circle is Kf/4. 
This property excludes exceptional configurations, appearing e.g. for the Wilson 
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action and causing serious problems in numerical simulations. As a consequence 
of eq. (Q) the solutions with ImA ^ come in complex conjugate pairs and satisfy 
(V'lTst/O = 0. Due to eqs. (|],||), those with A = are chiral eigenstates, 75^ = ±V- 
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